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Abstract 

In this paper we consider the three-particle density correlation function 
for a fractional quantum Hall liquid. The study of this object is motivated 
by recent experimental studies of fractional quantum Hall systems using 
inelastic light scattering and phonon absorption techniques. Symmetry 
properties of the correlation function are noted. An exact sum-rule is 
derived which this quantity must obey. This sum-rule is used to assess 
the convolution approximation that has been used to estimate the matrix 
elements for such experiments. PACS Numbers: 73.40.Hm, 73.20.Mf, 
72.10.Di 

1 Introduction 

The fractional quantum Hall effect has been one of the key topics in condensed 
matter physics for the last 15 years. While much of the experimental work on 
this phenomenon has been focussed on electrical transport measurements, re- 
cently attention has turned to the development of complementary spectroscopic 
techniques. In the case of both phonon]!], |2| and inelastic light scattering H || 
experiments the spectroscopic probe couples to the electronic density of the 
system and so naturally interacts with its collective density fluctuations. The 
theory of these collective modes was first developed by Girvin, MacDonald and 
Platzman (GMP)g in a manner analogous to that used by Feynman for the 
phonon and roton modes of superfluid helium. Unlike the case of liquid helium, 
there is no gapless phonon-like behaviour at small wave vector; the collective 
modes of the fractional quantum Hall systems are gapped at all wavelengths: a 
consequence of the incompressibility of these states. However, in common with 
liquid helium, the collective mode dispersion does have a well defined minimum, 
hence the quanta of collective excitations of the fractional quantum Hall states 
are referred to as magnetorotons. 

One can envision two types of process involved in both the phonon and 
light scattering experiments. Firstly a magnetoroton can be created when the 
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electron liquid absorbs one of the probe quanta (phonon or photon) , we refer to 
this as type A absorption. Secondly, an existing magnetoroton can be scattered 
into a different state by the absorption of a probe quantum, this we call a type B 
process. At zero temperature the electron liquid will be in its ground state (well 
described by Laughlin's wave function^) and there will be no magnetorotons 
present, hence only the type A process can occur. Because the magnetoroton 
dispersion is gapped at all wave vectors there is a threshold frequency for this 
process: no probe quantum can be absorbed whose energy is less than the 
minimum gap, A*. Type B processes will occur at any non-zero temperature 
and have no threshold: in principal, arbitrarily small energies can be transferred 
to the electron liquid by these processes. The theory of the Type A processes has 
been discussed by He and Platzman for the inelastic light-scattering case and 
by our group || for the phonon case. These type B processes will determine the 
leading finite-temperature corrections to the zero temperature theory outlined 
in ||. The theory of Type B processes is less well developed but potentially 
rather important for the understanding of the phonon experiments because these 
typically use pulses containing a black body distribution of phonon energies 
characterized by a pulse temperature < A* and use the temperature of the 
electron liquid itself to monitor the rate of phonon absorption. At an arbitrary 
electron temperature there will be competition between rare type A processes 
in which high energy phonons are absorbed and common type B processes in 
which low energy phonons are absorbed. In recent time- resolved experiments Q 
the change in electron temperature is recorded over the period of time in which 
the pulse is in contact with the electrons. At early times the electron liquid 
is at very low temperatures and so type A processes will dominate the energy 
transfer and hence the rate of change of electron temperature. At later times the 
electron temperature will be higher, because of all the earlier type A absorptions 
and type B processes will become comparable in determining the energy transfer 
rate. 

The purpose of this paper is to derive a sum rule obeyed by the matrix 
elements which control the rates for the type B processes and to use this to 
test the best existing approximation scheme for this object: the colvolution ap- 
proximation developed for studying the corresponding quantities for superfluid 
Helium-4 fl2|| . As discussed in these matrix elements are also important 
for the understanding of the effects of disorder on the form of the magnetoroton 
dispersion. In all of the following, natural units in which l c = yJh/eB — 1 and 
e 2 / AtteokIc = 1 (where k is the bulk dielectric constant of the material in which 
the 2des is formed) are used. 

2 Magnetorotons 

Girvin, MacDonald and Platzman|| supposed that the low-lying collective ex- 
citations of a fractional quantum Hall fluid could be obtained in a manner anal- 
ogous to that used by Feynman in his classic works on superfluid helium. They 
proposed that a collective excitation with 2d wave vector q would be described 
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by the wave function 



w- ( a) 

where |f2) is the ground state of the system (for which Laughlin has provided 
an excellent wave function, at least for the cases where the Landau level filling 
v is of the form 1 jm for m odd) and p q is an operator that is derived from the 
conventional density operator 

JV 



9q = 5>**« (2) 



i=l 

by projection onto the lowest Landau level. In the notation used by MacDonald 
in § 

N 

p q = X>(q) (3) 

i=l 

B (q) = Voe^Vo (4) 

where Vq acts on the Hilbert space of a single electron to project out the states 
within the lowest (spin-polarized) Landau level. Girvin and Jach|l(J investi- 
gated the properties of these projected operators and, in particular deduced 
that 

Bi (k) Bi (q) = e k *^ 2 i? 4 (k + q) (5) 

where k = k x + ik y . This leads, bearing in mind that the projection opera- 
tors for one particle commute with operators associated with the other, to the 
commutation relation for the projected density operators 

[p k ,P q ]-( e k * q/2 -^ k/2 )p k+q 

= i$(k,q)p k+q (6) 

where 

$(k,q) = -z( e k *^ 2 - e ^ k / 2 ) (7) 
= e kq / 2 2sinQkAq 

The normalization in equation [l] includes a matrix element which has the 
form of a projected static structure factor 
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GMP showed that this could be related to the true static structure factor 

1 

N 



s(q) = ^(n\ P ^ Pci \n) (9) 



via the simple relation 

s{q)-e-i 2 l 2 = s{q)-l . (10) 
The magnetoroton energy can be written in the form 

A(fc) = (k|«|k) - (n|w|n) 
_ (n\p_ k [n,p k ]\n) 



Ns(k) 

7(fc) 

s(k) 



(ii) 



and GMP showed that the projected oscillator strength, f (k), can be written 
explicitly as an integral involving s and $ only. Hence, obtaining s (k) from the 
Laughlin wave function, they found the dispersion relation for these excitations. 

Consider now a bosonic probe quantum (such as a phonon or photon) la- 
belled by a 3d wave vector Q which couples to the 2d electrons via a hamiltonian 
of the form 

flint = M Q<°q (°Q + a -q) ( 12 ) 
Q 

where <zq annihilates a probe quantum with the given wave vector, p q is the elec- 
tron density operator described above (un-projected) and Mq is some coupling 
function. For the detailed forms of Mq relevant to the experimental systems 
see j^, ^j. The probability per unit time for a magnetoroton initially in a state 
|k) to absorb a probe quantum with wave vector Q and be scattered into the 
state |k + q) is given by Fermi's golden rule as 

= y |M Q |2 1 (k + q| p i |k> I 2 6 (A (|k + q|) - A (fc) - (13) 

where wq is the energy of the probe quantum. This, as usual with lowest order 
perturbation theory, factorizes nicely into a part which depends on the well 
characterized details of the probe and the material in which the 2des is formed 
and a part that only involves the states of the correlated 2d electrons. The 
relevant matrix element in the latter is 

(k + q| /0q |k)= - P(k ' q) = (14) 
^s(\k + q\)s(k) 

where we have replaced /? q with its projected counterpart by virtue of the idem- 
potence of the projection operators and we define the three-point correlation 
function 

P(k,q) = l(r!|p_ k _ q p q p k |0) . (15) 
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3 Symmetry Properties of the Correlation Func- 
tion 

We can deduce two symmetry properties of P (k, q) straightforwardly. Firstly 
we can use the fact that p q = p_ q to deduce that 

{P(k,q)}* = |l(r!|p_ k _ q p q p k |r!) 

= ^(^|P-kP-qPk+q|^) 

= P(k + q,-q) . (16) 

Secondly we can use the commutation relation for the projected density 
operators, as derived by Girvin and Jach|l(|to simplify 

P (k, q) - P (q, k) = 1 (Q |p_ k _ q [p q , fcj | A) 

= (q, k) ^ (fi |p_ k -qPk+q| 

= t$(q 3 k)S(|k + q|) . (17) 

4 Derivation of the Sum Rules 

4.1 Structure Factor Sum- Rule 

In order to show the basic idea of the method we will begin with a simple case 
and prove a sum-rule for the projected structure factor itself. From the work of 
GMP H we know that 

s(q)-e-« 2 / 2 =s(q)-l . 
In first quantized notation we have that, for a system of TV particles 



in) 



Hence we have that 



Now 



J ( S (q)-l)d 2 q=l^(n| J e-^-^d\\Q) 
J e -*«-<?*-*i) < jP cl = 6(f i -¥ j ) 
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so that 



J ( S (q)-l)d 2 q=l^(fi|<5(r J -r J )|fi) . 

Inserting a representation of unity in terms of the complete set of iV-particle 
position eigenstates |{rj}) gives 

f ( s (q)-l)d 2 q=l]T f [Y[d 2 rA <Q| S (?,-¥ j) |{r0) ({rj |fi) 
J ijtj J \i=i ) 

= jjT,[ (n^riU^-r^^Kr.DUrOlO) 

i^j J \»=i / 

The only configurations that could contribute to this integral are those in which 
two particles co-incide. Such configurations, however, have zero weight because 
of the Pauli exclusion principle. Hence we conclude that 

|( S (q)-lKq = . 

From this we find, straightforwardly that 

J S (q) d 2 q - 2ir J™ e^^qdq = J (s (q) - 1) d 2 q = 



or 



J s (q) d 2 q = 2tt 

Of course, this result can be derived more simply by noting that 

s (q) = 1 + p j e-^ r (g (r) - 1) d 2 r + 4tt 2 P 5 2 (q) 

where g (r) is the radial distribution function. Hence 

J (s (q) -l)d 2 q = pj J e-^ r d 2 q (g (r) - 1) d 2 r + Att 2 P J 8 2 (q) d 2 q 

= 4tt 2 P J 6 2 (r) (g (r) - 1) d 2 r + in 2 : p 
= ATi 2 pg (0) 

The Pauli principle ensures that g (0) = so that we again recover our sum-rule. 
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4.2 Three Point Correlation Sum-Rule 



Now let us derive our principal result: a sum rule for the three point correlation 
function. As shown by MacDonald et al.p|, the three point correlation function 
can be written in the form 



P(k,q)=F (k,q) + ^ 3 ) (k,q) 



(18) 



where 



Po (k, q) = e -slql 2 -3k*q j e -^l 2 + s (| k |) _ X J 

+ e-^ k +^' k (s (|q|) - 1) + e^* k (s (|k + q|) - 1) 
Once again, k =k x + ik y is the complex representation of the vector k and 



^ (3) (k,q) = ^ E (« 



e -i(k+q) r; e »q r TO gik-r r . 



(19) 



is the un-projected three point correlation function of the quantum liquid. 
We wish to find a sum rule of the form 



i^P(k,q) = F(k) 



A 



(20) 



The term that is likely to cause us difficulty is the un-projected three point 
function, however 



l5>(3)(k,q) = l£l £ (0 



A^ N 



g -i(k+q)-ri e iq-r me ik-i 



n 



(21) 



The sum over q is simply a delta function, so that 

!]>>< 3 )(k,q) = l £ (n|e- ik - r ^ 2 ( ri -r m )e ik - r »|0) (22) 



but, as seen above, this must be identically zero, the only configurations for 
which the delta-function is non-zero are ones for which the ground state wave 
function vanishes by virtue of the Pauli exclusion principle. Hence ^ q (k, q) 
Vk. 

The remaining terms are simplified by noting that 



e -z k +s(k) -1 = s(k) 



(23) 
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and by removing possible singularities by writing 

s (q) - 1 = h (q) + NS^.o (24) 

where 

h(q) = p J(g(r)-l)e^d 2 v (25) 

and g (r) is the usual radial distribution function of the electron liquid. Hence 
we are left with 



I £ P (k, q) = 1 £ e-iW-i^ (|k|) + 2p e^l k l 2 
q q 



q q 

I^ e -^|q| 2 -^- q?( | k | ) + 2/ , oe -i|kP (26) 

q 

2e -il k l 3 l^e-^(|q|) . (27) 



A 
q 



As shown below in the appendix the summations can be carried out in the limit 
N, A — > oo, N/ A = p to give 



l^p(k, q )^iM + 2poe -^ + 2e -m> ( _ Po) (28 ) 

q 

s{k) 



(29) 



2tt 

which is our desired sum-rule. 

We can use the symmetry properties derived above to derive a second version 
of this. Consider 

^£( P (k, q) - P (q, k)) = »1 £ $ (q, k) 5 (|k + q|) 
q q 



^*(q-M)s(5) ■ (30) 



= l A 

q 

It is shown in the appendix that this integral vanishes so that 

iE*o,<o-*? . (31) 

k 

Finally the reality of the right-hand side of ^8] allows a further form of the 
sum-rule to be written as 

^ E P( k + q ,- q) = £M . (32) 
q 

The basic idea used here could be extended to consider sum-rules for higher 
order correlation functions were they to be of interest. 
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5 Assessment of the Convolution Approxima- 
tion 

In their work MacDonald et al. estimated the small q behaviour of P (k, q) 
by using the convolution approximation |fl2| for the 3-particle distribution func- 
tion 

n( 3 >(r,r',r")= £ (Q| 6 (r - r,) 6 (r' - r m ) 5 (r"-r n ) \Q) 
which can be written as 

n { p (r, r', r") = p 3 {1 + ft (r - r') + h (r' - r") + h (r" - r)} 

+ p 3 {h (r - r') h (r' - r") + ft (r' - r") ft (r" - r) + ft (r" - r) ft (r - /)} 

+ p 4 y h(r~R)h(r'-R)h (r" - R) d 2 K 
The three-point function, ft' 3 - 1 (k, q) , can be written as 

ft^ (k,q) = i y d 2 rd 2 r'd 2 r"e- l ( k+q ) r e lq r 'e 4k r "n^ (r,r',r") 

which gives, using the convolution approximation 

ftf (k,q) = N 2 5 k 

+ ft (*) ft (q) + [ft (*) + ft (q) + ft (fc) ft (g)] ft (|k + q|) 

This approximation improves on the standard Kirkwood decoupling in that it 
correctly gives the q — > limit as 

ftf) (k,0) = (iV-2)( S (fc)-l). 

Since Pq (k, q) saturates our sum-rule we expect that 

l^ft( 3 )(k,q)=0 . 
q 

Using the convolution form gives 

j 53 fti 3 ) (k, q) = pN8 k>0 + N5 k , j h (?) + 2 P h ( fc ) 
q q 

+ J E i h (*) ft (?) + ( ft (*) + h (?) + h ( fc ) ft («)) ft d k + 

q 

Now 

= -P 
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so that 



F W = T E ^ (k, q) = (1 + h (fc)) ±5>(ff)* d k + ql) 



q q 




This final integral requires a specific form for the pair correlation function, /i (g). 
For the primary fractional quantum Hall states (y = 1/m, for odd m) this is 
known from Monte-Carlo studies of the Laughlin wave function [^| which lead 
to the form 



where the c m co-efficients are tabulated in [|5|. We have estimated F (k) numer- 
ically for the case v — 1/3 and it is plotted in figure 1, along with a plot of 
s (k) /2ir for comparison. Clearly F (fc) is not identically zero and is not even 
negligible in comparison to s (fc) /2ir. Hence we deduce that this approximation 
is deficient. MacDonald et al.jLlJ themselves pointed out that this approxima- 
tion could well be unreliable as it does not correctly reflect the particle-hole 
symmetry (y — > 1 — v) of the system. 

6 Summary and Discussion 

In this short note we have derived a sum rule (with symmetry related variants) 
for the static structure factor and the three-point correlation function that will 
determine the leading finite-temperature corrections to the absorption rates of 
phonons and photons in spectroscopic studies of the fractional quantum Hall 
effect. In principle the approach followed would allow the formulation of similar 
sum-rules for higher order correlation functions should such objects ever be- 
come relevant to experimental work. The sum-rule for the three-point function 
has been used to assess the validity of the convolution approximation Jl2j for 
fractional quantum Hall systems and has shown that it does not capture all of 
the physics. 

This work was supported by the EPSRC (UK). 

A Evaluation of Integrals 

First of all we need to evaluate the following integral for later use: 



h(q) 



= —ve 



in 
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now ll^ 



rZ7T 

/ e" 2COS *cos(zsin ( /.)d0 = 27r 
Jo 

p27T 

/ e- zcos ^sin(2siii0)# = O 
Jo 



so that 



2tt 

e- ze **# = 27r Vz . (34) 



Now we need to evaluate 



in 



the limit .4 — > oo, TV — > oo, iV/„4 = p this becomes 



1 /*°° /•27T -i 

-L/ ^e-W #e-l^ = i- . (36) 



4tt 2 7 y 2^ 



Similarly 



1 f 00 

= — / dqqh(q) . (37) 

The final integral can be evaluated quite generally. 

J dqqh 0) = ^ J d 2 qh (q) 

= 2irp Q h (0) . (38) 

Now for any spin-polarized fermi system g (r) = 1 + h (r) -^0asr->0asa 
consequence of the Pauli principle, hence h(0) = — 1 and 

1 f°° 

— J dqqh (q) = -p . (39) 
Finally we need to evaluate 

i i r°° /" 2Tr / 1 

-^(q-k.kjafo)-^ e " fe2/2 y ^^7 ^e^^sinUfcgsm, 

(40) 

As we have already seen the 4> integral vanishes and so, therefore does the whole 
expression. 
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B Figure Caption 



Figure 1: A plot of the function F (k) estimated numerically (full line) with a 
plot of s(k) /2tt plotted (dashed line) for comparison. 
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